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Abstract. Planetary nebulae (PNe) have the potential to revolutionize our understanding of 
extragalactic stellar populations. Indeed, in many systems, bright PNe are the only individual 
] objects identifiable from the ground, and, even more often, they are the only stars that are 
. amenable to spectroscopy. We review the ways in which ensembles of PNe can be used to probe 

■ the metallicity, age, and history of a stellar population. In particular, we discuss three cases: the 
weak line spectroscopic regime, where one has knowledge of the line-strengths of faint forbidden 
lines such as [O III] A4363, a bright-line regime, where only the strongest emission lines are 
visible, and the photometric regime, where the only information available is the [O III] A5007 

<^ , luminosity function. We show that each of these cases, when properly calibrated, can provide 

■ unique insights into the objects that make up a stellar population. 
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1. Introduction 

■ In order to understand how galaxies are put together, we must be able to characterize 
the collective properties of their stars. Information about a system's age (or history of 
star formation), metaUicity (or history of chemical enrichment), initial mass function, 

O-i and binarity are aU critical for modehng the processes which define galactic evolution. 
Q Unfortunately, due to the great distances involved, obtaining such data is often extremely 
^ difficult. 

c/2 ■ This is where planetary nebula (PN) research can help. All stellar populations older 
[ than ~ 10® yr make planetary nebulae, and, in a large galaxy, hundreds of PNe are 

!>■ available for study. More importantly, PNe are bright. In systems outside the Local 
k>( Group, PNe are often the only individual stars that can be identified from the ground, and 

5_j are almost always the only (non exploding) stars that can be studied spectroscopically. 

Ci ^ Thus, planetary nebulae can provide unique insights into the history of a stellar system. 

The amount of information provided by a population of PNe depends on the type of 
data that one can acquire. In nearby systems, it is possible to use 4-m and 8-m class 
telescopes to observe all the important emission lines produced by a PN, including the 
critical diagnostic lines of [O III] A4363, [S II] AA6716,6731, and [N II] A6548,6584. 
With this, one can directly measure a population's abundance distribution function and 
constrain its history of chemical enrichment. In more distant galaxies, these weak emission 
features are usually unavailable, but one can still use brighter lines such as [O III] A5007, 
[Ne III] A3869, Ha, and H/3 to probe stellar content. Finally, if PN spectroscopy is 
impossible, one can still investigate galactic history via the [O III] A5007 planetary nebula 
luminosity function (PNLF) . We discuss each of these regimes below. 
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2. The Weak Line Regime 

A decade ago, measurements of the critically important [O III] A4363 emission line 
were mostly restricted to the PNe of the Galaxy and the Magellanic Clouds. For ex- 
ample, it was only through heroic efforts that Richer, Stasihska, & McCall (1999) and 
Jacoby fc Ciardullo (1999] were able to detect this feature in a sizeable sample (42) of 



PNe in M31's bulge, and thereby measure the population's oxygen abundance. These 
authors obtained the interesting result that the metallicity spread for M31's bulge stars 
is over a dex, and that the region's mean oxygen abundance is slightly less than solar. 
Unfortunately, given the limitations of the 4-m apertures that were available at the time, 
follow-up studies proved to be impossible. 

Today in this era of 8-m class telescopes, it is possible to obtain abundances for hun- 
dreds of PNe throughout the largest of the Local Group galaxies, M31 and M33. With 
these data, one can characterize the systems' abundance distribution functions, not only 
for oxygen, but for the a-process elements of neon, sulfur and argon, all of which are 
unaffected by stellar evolution. Moreover, because one knows the distances to these ex- 
tragalactic PNe, the combination of [O III] A5007 photometry and spectrophotometry 
can be used to derive the objects' absolute line fluxes. This is important, because from 
this information, it is possible to deduce the locations of the central stars on the HR 
diagram, estimate their masses (via comparisons with post-AGB evolutionary tracks; 
Vassiliadis & Wood 1994), and infer the masses of the stars' progenitors (via the initial 
mass-final mass relation; Weidemann 2000). In other words, deep spectrophotometry of 
a large number of planetary nebulae should allow us to trace the chemical enrichment of 
a stellar population back through time. 

Although such an analysis is exciting, it is also extremely difficult. For example, con- 
sider the PNe of the Large Magellanic Cloud. Excellent spectrophotometry and nebular 
models exist for 60 objects from the work of Dopita et al. (1991a, b; 1997). If one 
uses these data to derive a-element abundances, and central star luminosities, tempera- 
tures, and (ultimately) masses, then a clear trend emerges: PNe from lower-mass (older) 
progenitors are systematically more metal-poor than their younger counterparts. This 
suggests that we can, indeed, use PNe to probe a galaxy's history of chemical enrich- 
ment. However, if one instead uses ultraviolet measurements of the nebular and stellar 
continuum to place the central stars on the HR diagram (Herald & Bianchi 2004; Villaver 
et al. 2003; Stasihska et al. 2004), then much of this correlation goes away. Clearly, in 
order to exploit PNe to their fullest, our inferences about central star properties need to 
be made more robust. 



3. The Bright-Line Regime 

Once outside the Local Group, the [O HI] A4363 emission line of planetary nebulae 
becomes extremely hard, if not impossible, to detect. Although this makes traditional 
nebular modeling impossible, it does not preclude the use of PNe for population studies. 
Indeed, ,Dopita, Jacoby, fc Vassiliadis (1992), and Mendez et al. (2005) have shown that 



the ratio of [O HI] A5007 to H/3 can be used to place useful constraints on a PN's oxygen 
abundance. 

But one can do better than that. Due to the poor distances of Galactic PNe, samples 
of Milky Way planetaries are extremely heterogeneous in nature. With high and low-core 
mass PNe mixed together, it is nearly impossible to separate out planetaries at their peak 
luminosity from PNe that have already faded by several magnitudes. In extragalactic 
systems, this is not a problem. Since the absolute luminosities of all the PNe are known. 
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Figure 1. The PNc of the LMC (left) and the large elliptical NGC 4697 (right) plotted in a 
three-dimensional phase space, with the absolute [O III] magnitudes displayed on the y-axis. 
Note how this type of plot clearly distinguishes the different PN populations. 

it is possible to define statistically complete samples of objects at the bright end of 
the luminosity function and be assured that all are in a similar phase of evolution. For 
example, PNe in the Milky Way exhibit a wide range of excitation, with [O III] A5007 
to Ha emission line ratios (R) ranging from ^ 4 to 0.1. However, when one studies 
extragalactic PNe, a pattern emerges: the brightest PNe all have R > 2, regardless of 
their parent stellar population (Ciardullo et al. 2002). 

In fact, the best way to analyze the bright emission lines from a homogeneous PN 
dataset is through a differential, multi-dimensional analysis. A schematic example of this 
approach is shown in Figure 1. In the figure, PNe of the LMC (from Meatheringham 
& Dopita 1991a, b) and NGC 4697 (Mcndez et al. 2005) are plotted in a 3-dimensional 
space that has the absolute [O III] A5007 line flux on the y-axis, and the [O III]/H/3, and 
[Ne III] A3869/H/3 line ratios on the x- and ^;-axes. The differences in the populations are 
obvious. With a reasonable amount of effort, this type of data could be acquired for a 
range of stellar populations out to ~ 10 Mpc, and serve as the basis for an investigation 
of PN production and evolution as a function of galactic environment. 



4. The Photometric Regime 

Spectroscopy of extragalactic PNe is time consuming, and in many cases impossible. 
Hence one often must rely on [O III] A5007 (or possibly Hq) photometry to probe a 
PN population. Observations in over 50 galaxies have shown that the bright-end of the 
[O III] A5007 planetary nebula luminosity function is well-described by the truncated 
exponential 

7V(M)oce°-307M|i_e3(M-M*)| 

where 

M = -2.5 log F5007 - 13.74 (4.2) 

and M* = —4.47 is the absolute magnitude of the most luminous planetary (Jacoby 
et al. 1992; Ciardullo et al. 2002). Since M* is an excellent standard candle for all but 
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Figure 2. The [O III] A5007 planetary nebula luminosity functions for M31's bulge, M33's disk, 
and the Small Magellanic Cloud. Note the differences between the populations: star forming 
systems have a strong deficit of intermediate-brightness planetaries. 



the most metal-poor systems, its brightness cannot be used for stellar population studies. 
However, two properties of the PNLF do hold promise for population measurements. 

The first of these is a "dip" in the luminosity function which starts ^ 2 mag below 
M*. This feature is not present in every galaxy: in M31, the faint-end of the PNLF is 



well described by the simple ,Henize &: Westerlund (1963) exponential, which is expected 
from an ensemble of uniformly expanding spheres ionized by non-evolving central stars 
(CiarduUo et al. 1989, 2002). However, as Figure 2 illustrates, star-forming systems such 
as the Small Magellanic Cloud and M33, have a factor of ~ 2 fewer PNe than predicted 
from the exponential law (Jacoby & De Marco 2002; CiarduUo et al. 2004). Clearly, this 
part of the PNLF is sensitive to some aspect of stellar population. 

The second population dependent feature of the [O III] PN luminosity function is its 
normalization. Following CiarduUo et al. (1989)[ this property is usually parameterized 
by a, the number of PNe detected above a certain absolute magnitude (typically 0.5 or 
2.5 mag below M*), divided by the total (bolometric) luminosity of the stellar parent 
population. In the star-forming galaxies of the Local Group, a is in good agreement with 
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the predictions of stellar evolution models (Buzzoni, Arnaboldi, & Corradi 2006), but in 
older stellar populations, this is not always the case. In fact, some elliptical galaxies have 
almost an order of magnitude fewer PNe than one might expect. This variation has the 
potential to be an extremely powerful probe of early-type systems. 



5. Explaining the PNLF 

To date, there has been no attempt to simultaneously model both the shape and 
normalization of the PNLF. Moreover, those models which have concentrated solely on 
the shape of the luminosity function have failed spectacularly, either by requiring a 
constant star-formation rate for elliptical galaxies (Mendez & Soffner 1997) or predicting 
a bright-end cutoff that is several magnitudes fainter than what is observed (Marigo et al. 
2004). Obviously, an important piece of physics is missing from these analyses. 

The main difficulty in modeling the PNLF comes from the absolute magnitude of its 
cutoff. An M* PNe emits ^ GOOL© in its [O III] A5007 emission line. Since no more than 
~ 10% of the central star's ffux can come out in this line, the luminosity of the source 
powering the nebula must be at least ~ 6, OOOL©. Unfortunately, the only post-AGB cores 
capable of generating this much energy are those with masses of M > O.BM© (Vassiliadis 
& Wood 1994), and these require progenitor stars with M > 2Mq (Weidemann 2000). 
Such stars have lifetimes of less than a Gyr, and are not present in any number in elliptical 
galaxies {e.g., Trager et al. 2005). 

As Figure 3 demonstrates, the situation is even worse than it looks. Spectroscopic 
observations in the Large Magellanic Cloud (Meatheringham & Dopita 1991a, b) and M31 
(Jacoby & CiarduUo 1999) demonstrate that internal extinction is extremely important 
in [O III]-bright planetary nebulae. While no PN is observed to have an [O III] luminosity 
brighter than Af*, many PNe have intrinsic [O III] fluxes that are much more luminous 
than this, sometimes by a factor of ^ 2. Such a high luminosity can only be generated 
by a massive post-AGB core. 

In star-forming systems, the presence of high core-mass central stars is not a problem. 
But what produces these objects in old stellar populations? The number of PNe observed 
in a galaxy, normalized to the galaxy's underlying (bolometric) luminosity is given by 

a^B-T-f (5.1) 

where B is the population's luminosity specific stellar evolution flux, r, the lifetime of the 
PN stage, and /, the fraction of stars which evolve through the PN phase. The value of B 
is known from the theory of stellar energy generation: all old stellar populations, regard- 
less of their age, metallicity, or initial mass function, have B 2 x 10^^^ stars yr~^ Lg^ 
(Renzini & Buzzoni 1986). The lifetime, r, is similarly well-constrained, as models show 
that [O III]-bright PNe spend ~ 500 yr in the top 0.5 mag of the luminosity function 
(Marigo et al. 2004). Therefore, observed values of a in the range ^ 5 x 10~^° < ao.5 < 
3 X 10~^ PNe Lq^ (CiarduUo et al. 2005) translate into stellar fractions of between 5% 
and 30% (see Figure 4). Since absorption line spectroscopy of early- type galaxies place 
the fraction of young stars far below this value (Trager et al. 2005), the bright PNe of 
these systems cannot be evolving via normal (single star) stellar evolution. We must look 
elsewhere for their explanation. 

As of this date, two scenarios have been proposed to explain the existence of [O III]- 
bright PNe in early- type galaxies. The first, described by CiarduUo et al. (2005)[ involves 



binary coalescence and blue stragglers. Although elliptical galaxies possess few stars with 
M ~ 2Mq, they do contain many stars with M ^ IMq. If some of these stars come 
into contact while still in their hydrogen burning phase, the result may be conservative 
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Figure 3. The observed ratio of [O III] A5007 to Ha plotted against [O III] A5007 absolute mag- 
nitude for those bright PNe in the LMC and M31 with published spectrophotometry. The arrows 
display the effects of internal extinction, as estimated from the objects' Balmer decrement. The 
location of the PNLF cutoff, M*, is identified via the dotted line. Note that, although no FN 
is ever observed to have an [O III] A5007 magnitude brighter than M*, many have intrinsic 
luminosities much brighter than this. 



McCrea (1964]} style mass transfer and coalescence. This mechanism, which has been 
associated with the creation of blue stragglers (Mateo et al. 1990; Carney et al. 2001), has 
the potential to build the requisite high mass cores, even in systems as old as 10 Gyr. 

The blue straggler hypothesis has several advantages. The first, and most obvious, is 
that such stars are known to exist. In fact, if blue stragglers live ~ 5 x 10® yr (Lombardi 
et al. 2002), then the ratio of [O III]-bright PNe to blue stragglers in systems such as 
the Ursa Minor dwarf galaxy and the Milky Way halo is approximately the same as that 
of the objects' lifetimes (CiarduUo et al. 2005). This is consistent with the idea that one 
object evolves into the other. Moreover, the binary merger hypothesis provides a natural 
explanation for the observed inverse correlation between a and the excess ultraviolet 
flux often associated with elliptical galaxies. As a stellar system ages, its main sequence 
turnoff mass decreases, and this decreases the mass of its post-AGB cores (Weidemann). 
In the extreme, these low-mass cores will evolve too slowly to form planetary nebulae, 
and instead produce naked UV-bright sources (Greggio & Renzini 2000). Meanwhile, the 
smaller turnoff mass also reduces the number systems capable of merging into stars above 
the critical ~ 2Mq threshold. Consequently the number of M* planetaries will also drop. 
If this scenario is correct, then a can be used in conjunction with a galaxy's UV upturn 
to probe the precise location of a population's main-sequence turnoff, even when this 
turnoff is far below the limit of detectability. 

The blue straggler hypothesis only runs into problems when one considers the con- 
straints it imposes on the integrated spectra of elliptical galaxies. As pointed out by 
Xin fc Deng (2005)| a sizeable population of blue stragglers can mimic the presence of 
young stars, and lead to a significant underestimate for the age of the stellar population. 
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Figure 4. The bolometric-luminosity specific number of PNe in tlie top 0.5 mag of tlie [O III] 
PNLF plotted against tiie parent galaxy's optical and ultraviolet color. Circles indicate elliptical 
galaxies, squares are lenticulars, triangles are spiral bulges, and crosses represent systems which 
contain obvious evidence of star formation. The error bars represent formal 1 a uncertainties 
and are lower limits to the true errors. The right hand axis translates the FN densities to the 
fraction of stars turning off the main sequence which evolve into [O III]-bright planetaries, under 
the assumption that the PNe are products of single star evolution. 



In fact, Buzzoni, Arnaboldi, & Corradi (2006) suggest that any blue straggler contribu- 



tion to an old stellar system must be considerably less than 10%, in order to avoid 
measurable effects on the galaxy's colors. However, whether this is a serious problem for 
models of elliptical galaxies is still an open question. 

An alternative scenario for the presence of M* [O III] A5007 sources in early-type 



populations involves symbiotic stars. Soker(2006) has proposed that mass transfer from 



a proto-planetary post-AGB star onto a companion white dwarf can produce the energy 
and emission lines needed to mimic a bright planetary nebula. This model, which has 
many commonalities with models of super-soft x-ray sources (albeit with a much higher 
accretion rate), has the advantage that it does not stress population synthesis predictions: 
since the lifetime of this accretion phase can be much longer than ~ 500 yr, far fewer 
[O III] systems are needed to populate the bright end of the luminosity function. However, 
the model has no natural explanation for the similarity between the maximum brightness 
of a true planetary nebula and one of these symbiotic sources, nor does it explain the 
inverse correlation between a and a population's UV upturn. 

If there is indeed a second population of objects that evolve through a non-traditional 
(binary-star) mode of stellar evolution, then the dip in the PNLF has a simple explana- 
tion. High core mass planetaries, whether from single star or binary evolution, will always 
form a luminosity function that has a distinctive dip at intermediate magnitudes. The 
only way to avoid such a feature is to fill in the dip with lower mass objects. The ampli- 
tude of the dip, along with its precise location, would then be a measure of how many of 
these difhcult-to-observe low-mass stars are present in the population and contributing 
to the PN luminosity function. If this idea is correct, then the agreement between M31's 



PNLF and the Henize & Westerlund (1963) exponential law is purely coincidental 
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